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THE CONTACT PROCESS ON THE COMPLETE GRAPH WITH 
RANDOM VERTEX-DEPENDENT INFECTION RATES 

JONATHON PETERSON 



Abstract. We study the contact process on the complete graph on n vertices where the rate 

at which the infection travels along the edge connecting vertices i and j is equal to XwiWj/n 

for some A > 0, where Wi are i.i.d. vertex weights. We show that when -B[it)i] < cxo there is 

a phase transition at Ac > so that for A < Ac the contact process dies out in logarithmic 

time, and for A > Ac the contact process lives for an exponential amount of time. Moreover, 

we give a formula for Ac and when A > Ac we are able to give precise approximations for the 

probability a given vertex is infected in the quasi-stationary distribution. 

Q^ ' Our results are consistent with a non-rigorous mean-field analysis of the model. This is in 

Q, , contrast to some recent results for the contact process on power law random graphs where the 

• ' mean-field calculations suggested that Ac > when in fact Ac = 0. 



1. Introduction 

The contact process is a simple model for the spread of a disease. The standard model of 
the contact process on a graph G = {V, E) is described informally as follows. Fix a parameter 
A > and a set of vertices A dV . At time t = only the vertices in A are infected. As time 

op ■ progresses, each uninfected vertex x becomes infected at rate equal to A times the number of 

currently infected neighbors, and each infected vertex becomes a healthy (uninfected) vertex 

jn I at rate 1. More formally, the contact process is a continuous time Markov process rn € {0, 1} 

^^ ■ with generator 



>C/(r?) = Y. ^(^) + (1 - ^(^))^ E '?(y) (/('?") - /('?)) ' 
ij x&V \ y^x / 

^ ■ where / : {0, 1} — )• M is a bounded function, x ^ y means that vertices x and y are connected 

by an edge in E, and rj^ is the configuration obtained from t] by switching the value of r]{x). 
That is, 

.'(.) = If ,, ^^" (1) 

I 1 — r][z) z = X. 

The contact process is also sometimes referred to as the susceptible-infected-susceptible (SIS) 
epidemic model. 

The behavior of the contact process depends on the parameter A, and as A increases the 
infection spreads faster and it takes a longer amount of time for the contact process to die 
out (i.e., reach the absorbing state of all healthy vertices). It is then natural to ask if there 
is a critical values of A for which the contact process exhibits a phase transition. The contact 
process on the integer lattice Z has been well studied, and it is known that there is a Ac > 
such that the contact process started with a single vertex infected dies out with probability one 
if A < Ac and survives forever with positive probability if A > Ac |Lig99| . On any finite graph the 



Date: May 3, 2010. 

2000 Mathematics Subject Classification. Primary: 60K35; Secondary: 60K37, 82B26, 05C80. 

Key words and phrases. Contact process, random environment, phase transition. 

J. Peterson was partially supported by National Science Foundation grant DMS-0802942. 
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contact process always eventually dies out, and thus it is not immediately clear how to define a 
phase transition. However, for the contact process on [— ?i, n] C Z it is known that for A < Ac 
(where Ac is the critical parameter for the contact process on Z ) the contact process dies out 
by time C log n with high probability, whereas for A > Ac the contact process survives for time 
expjcn } with positive probability |Lig99 . In general, one says that the contact process on a 



family of finite graphs is sub-critical if the time until the infection dies out is logarithmic in 
the number of vertices and is super-critical if with positive probability the infection survives 
for a time that is larger than any polynomial in the number of vertices of the graph. The 
critical value Ac then identifies the phase transition of the contact process from sub-critical to 
super-critical. 

Recently, Chatterjee and Durrett |CD09] and Berger, Borgs, Chayes and Saberi |BBCS05] 
considered the contact process on two different models of power-law random graphs. A power- 
law random graph is a general term denoting a class of graphs where the distribution of the 
degree of a typical vertex has tails that decay like Cfc~"~^ as fc — )■ oo for some a > 1 and 
C > 0. Physicists had previously studied the contact process on power-law random graphs 
using non-rigorous mean-field calculations and concluded that if a > 3 then there was a critical 
value Ac > identifying a phase transition [PSVOlbl PSVOla]. However, for the two types of 
power law random graphs studied in |BBCS05] and [CD09j it was shown that in fact Ac = 
(i.e., the contact process survives for a long time for any A > 0). The long time survival 
of the contact process implies the existence of a quasi-stationary distribution. The mean- 
field calculations suggest that the average density p(A) of infected sites in the quasi-stationary 
distribution satisfies p{X) ~ CX^ for some /3 > as A — )• A+. However, upper and lower bounds 
on p{\) calculated in |CD09] show that the "critical exponent" /3 must be different from the 
mean-field predictions. 

In this paper, instead of studying the contact process on a random graph, we use a deter- 
ministic graph (the complete graph il'„) and instead make the rates at which the infections 
travel along edges to be random. We will choose the random infection rates in a manner that 
is inspired by the power-law random graph model of Chung and Lu |CL03j . It turns out that 
in this model the mean field calculations actually turn out to be correct. We are able to obtain 
a formula for Ac and approximations for the average density p{X) of infected sites that agree 
with the mean- field predictions. 

1.1. Description of the Model. We now turn to a description of the specific process studied 
in this paper. We will denote the set of vertices of Kn by [n] := {1, 2, . . . n}. Let w = {wi}'^-^ 
be a sequence of non-negative real numbers, and assign weights w„ = {wi}^^^ to the vertices 
of Kn- Then, given the vertex weights, the rate at which infections are transmitted from i to 
j (or j to i) is equal to XwiWj/n. The rate at which infected vertices become healthy is kept 
constant at 1. Formally, for fixed w„, A > 0, and A C [n] we let r]t be the Markov process with 
r/o(-) = Ia(-) and generator 

(^w„/)(r/) = f; ( r?(.) + (1 - r?(i)) J] r?(j)^^ ) Uiv') " fiv)), 

where rj^ is defined as in ([1]) . The law of rft in this case will be denoted by PZ' ■ Corresponding 
expectations will be denoted by E^ . Often we will be interested in the cases where either all 
vertices are initially infected or just a single vertex is initially infected. Thus, for notational 
convenience we will abbreviate Pw and P^ by P^ and Pw'*, respectively. Corresponding 
expectations will be denoted similarly. 

We will make the following assumption on the vertex weights 
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Assumption 1. The sequence of vertex weights w = {wi}^^ is an i.i.d. sequence of random 
variables with common distribution /i. 

Expectations with respect to the measure // on the vertex weights wiU be denoted E^. 

Remark 1.1. The essential property that we use is that hm.„_j.oo ;^ X^"=i /(u^O = E^[f{wi)\ 
for certain rational functions f{w). Thus, the statements of the theorems and the proofs 
remain essentially unchanged by letting w be an ergodic sequence or even a "nice" deterministic 
sequence (c.f. [MR98J ) with empirical distributions approximating the measure /u. 

Before stating our main results, we explain briefly the motivation behind this model. One way 
of constructing the contact process t/j with random edge weights described above is to construct 
for each i a Poisson point process Ni{t) with rate 1 and for each pair i ^ j a Poisson point 
process Nij{t) with rate XwiWj/n, such that all the Poisson point processes are independent. 
At each jump time of the process Aj(-), if the vertex i is infected it becomes healthy, and at 
each jump time of the process Ajj(-) if exactly one of the vertices i or j is infected then the 
other vertex becomes infected as well. 

If we consider the vertices to represent computers or individuals in a network then the jump 
times of the process Ajj(-) represent connections made between the respective individuals or 
computers (such as human contact or e-mail message sent). Thus, if we only keep track of 
connections formed over a short time period [0,(5], the resulting random graph will have the 
edge between i and j present with probability 1 — exp{—6XwiWj/n} w 6XwiWj/n. Chung and 
Lu |CL03j studied a model for random graphs where, given a sequence of vertex weights Wi, the 
probability that there is an edge connecting vertices i and j is proportional to WiWj/n. If the 
vertex weights Wi have power-law tails, then the degree distribution of the resulting random 
graph has power law tails. Thus, if the distribution fi on vertex weights has power law tails we 
might expect the contact process on Kn described above to be similar to the contact process on 
a power law random graph. It is therefore somewhat surprising that the behavior is different 
from the contact process studied on the power-law random graph models in |CD09l IBBCS05] 
and that the mean- field predictions in this case are actually correct. 

1.2. Main Results. We now turn to the statements of the main results of the paper. We 
introduce these via the mean-field heuristics which help to explain them. Assume for now that 
there are only a finite number of weights and that the number of vertices of weight x is exactly 
//({x})n. Let pt{{x}) be the probability that a vertex of weight x is infected at time t, and let 
Nt{x) = pt{x) iJ,{{x})n be the number of vertices of weight x that are infected at time t. Then, 
under the mean field assumption that ??i(i) and rjt{j) are independent (which is not true) we 
obtain that 

^^ = -N,{x) + Y.{^{{x})n - N,{x))N,{y)^. 
y 
Recalling that pt{x) = Nt{x) / [ii[{x})n) we obtain that 

dpt{x) 



-pt{x) + ^{l - pt{x))pt{y)Xxy fi{{y}) 



pt{x) + / (1 - pt{x))pt{y)Xxy fi{dy). (2) 



dt 



The above mean-field equation should also hold when the distribution of vertex weights is 
continuous as well. 

If the contact process survives for a long time then there should be a quasi-stationary distri- 
bution. Thus, we look for a stationary solution to dJD- That is, we want to find a function p(x) 
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such that 

p{x) = / (1 - p{x))piy)Xxy fi{dy) = / p{y)Xxy fi{dy) - p{x) / p{y)\xy ^i{dy). (3) 

Solving for p{x) we obtain that 

(TX I 

P{^) = T"^ > where a = X yp{y)fi{dy). (4) 

l + ax J 

Substituting the formula for p{x) on the left into the equation on the right yields the equation 

We wish to characterize for what values of A (depending on the distribution /i of vertex weights) 
there is a o" > which solves ([5]). Note that the right hand side of ([5]) is decreasing in a. Thus 
if there is a solution to ([5]) it is unique. Moreover, (assuming E^wi < oo) 

,2 



and 



lim A / u(dv) = 

^^^ J l + ay^"- ^' 



0-^0+ J 1 + cry 
This leads us to the following definitions. Let 



li™ / ^ 1 , -.. /^(^^) = / ^y"^ f^idy) = \E^{wl). 



{-ET-^ if EnvJi < oo 
e;^i ^ 1 (6) 

otherwise, 



and for A > Ac define a{X) as follows. 



(t(A) is the unique o" > that solves 1 = XEi 



M 



Wi 



1 + awi 



(7) 



Our main results are a confirmation of the above mean- field heuristics. The first result verifies 
the existence of a phase transition at Ac. 

Theorem 1.2. Let the vertex weights w be i.i.d. with distribution /i, and let Ac be defined as 
in dSD. 
(i) If X < Ac, then there exists a constant C > such that 

lim P^(r/ciogn = 0) = 1, fi-a.s. 

n— s>oo 

(ii) If X > Ac, then there exists a constant c > such that 

lim P^(r?ec" / 0) = 1, fi- a.s. 

n— >oo 

Our second main result states that when the contact process is super-critical, the probability 
of a vertex being infected may be uniformly approximated by the mean-field prediction ([3]). 

Theorem 1.3. Let the vertex weights w be i.i.d. with distribution fx, and let Ac and o"(A) be 
defined as in (^ and (UD, respectively. Then, for any A > Ac and e > there exist constants 
C, c > (depending on e) so that 

a{X)wi 



lim sup sup sup 

n-s^oo ie[Clogn,e':"] «e[n.] 



P^imii) = 1) 



1 + a{X)wi 



< e, n — a.s. 



Remark 1.4. We are actually able to prove an upper bound on the probability of a vertex being 
infected that is slightly better than what is implied by the statement of Theorem 11.31 See 
Proposition 14.31 for a precise statement. 
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In their study of the contact process on a random graph, Chatterjee and Durrett |CD09] 
analyzed the contact process at time t = e^"". The time e^" is large enough for the contact 
process to have stabilized but small enough so that with high probability it has not died out. 
Thus, the distribution of the contact process at time t = e^ when started with all vertices 
infected was called the quasi-stationary distribution. As in |CD09j . we define PnW to be the 
expected number of infected vertices at time e^ when initially all vertices are infected. That 
is, 



/j„(A) = E^ 



1 " 



v^{^) 



11 



«=i 



Then, we obtain the following simple Corollary of Theorem 11.31 

Corollary 1.5. If X > Ac, then 

a{X)wi 



/3(A) := lim /9„(A) = Ef, 



l + o-(A)u;i 



As mentioned above, upper and lower bounds on the critical exponent of Pn(A) were derived 
in |CD09j . Since 

E^ -— ^ E^[wi]a, aso--^0+, 

[1 + awij 

Cor ollarv 1 1 . 5 1 implies that p{X) and o"(A) have the same critical exponent. The following Propo- 
sition allows us to not only identify the critical exponent of o"(A) but also the leading constants. 

Proposition 1.6. Let p{wi > x) ~ Cx~^'^~^' as x — )■ oo for some C > and a > 2, and let 
Xc and (t(A) be defined as in ([6]) and ([7]). Then, as 5 ^ 0~^ , 

l/{3-a) 



a{6) 

log a{6) 

(t{Xc + 6) 

(t{Xc + 6) 

(j{Xc + 6) 



r CTTJa-l) 
\ sin(7ra) 
-1 
2C6 

/ — sin(7ra) 
VCA2(a-l)7r 

1 6 

3CA2 log(l/5) 
1 



^l/{3-a) 



l/{a-3) 



jl/{a-3) 



xm 



[wf] 



if a e (2,3) 
if a = 3 

*/aG(3,4) 

if a = A 

if a> A. 



Moreover, the conclusion in the case a > A also holds under the assumption that E^lwf] < oo 
(i.e., without any assumption on tail asymptotics of p). 

Remark 1.7. Since Theorem 11.31 shows that the infection probabilities agree with the mean- 
field predictions, these critical exponents agree with the mean-field predictions in [PSVOlaj . 
However, the case a = 4 was not considered in [PSVOlaj . 

Remark 1.8. Proposition [L6] states that in the case a = 3, a{d) = exp{ — ^(2^ + 0(1))}. Riordan 
|Rio05j has previously observed a similar fast rate of decay for the fraction of vertices in the 
giant component of percolation on the Barbasi- Albert model of power law random graphs (which 
have a = 3). 

Before proceeding to the proofs of the above results, we recall two useful facts about the 
contact process that we will use. Both of these facts are found in |Lig99| and can be proved 
using the what is called the "graphical representation." 
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• Monotonicity in infection rates. Increasing the infection rate along any edge only 
increases the number of vertices infected at any time t. In particular, if w < w in the 
sense that Wi < Wi for all i, an easy coupling argument implies that 

w < w ^ P^'^(7?t > 1b) < PZ'^ivt > 1b), V^, Bc[n\. (8) 

• Self-duality. For any subsets of vertices A and B, the probability that an element of 
B is infected at time t when the infection starts from A is equal to the probability that 
an element of B is infected at time t when the infection starts from B. That is, 

P^'^{i]t{i) = 1, for some i £ B) = P^'^{rit{i) = 1, for some i G A). 

An important special case of this is when A = {i} and B = [n]. In this case, we obtain 
that 

PZ''{m^o) = PZ{m{i) = i)- (9) 

The remainder of the paper is organized as follows. In Section [2] we consider the special case 
when the distribution fi has finite support. In this case we identify the contact process with 
a multi-dimensional birth-death chain. By analyzing the birth-death chain we then show that 
when A > Ac the contact process survives for time e*^" with high probability, and we also obtain 
a lower bound on the probability a given vertex is infected at any time t < e*^"". In Section [3] 
we extend the results of Section [2] to the general case by approximating // by a measure with 
finite support. Next in Section [4] we complete the proofs of Theorems 11.21 and 11.31 by proving 
the complementary results to those obtained in Sections [2] and [3l That is, we show that the 
contact process dies out with high probability by time Clogn when A < Ac, and we obtain 
upper bounds on the probability a given vertex is infected at any time t > C log n when A > Ac. 
The main technique used in Section [5] is a comparison of the contact process with a related 
multi-type branching process. Finally, the proof of Proposition 11.61 is given in Section [5j 

Acknowledgement. Many thanks to Rick Durrett for the many helpful discussions in the 
process of writing this paper. 

2. Finitely many types 

In this section we will analyze the above model under the assumption that the support of 
the /i on vertex weights Wi is finite. Let W = {Wi, W2, ■ ■ ■ , Wm) be the possible vertex weights, 
each occuring with probability pi = ^{wi = Wi). We will classify the vertices according to the 
weight that they are assigned. That is, a vertex with weight Wi will be referred to as a type 
Wi vertex. Since all type Wi vertices are equivalent, we need only to keep track of the number 
Xi[t) of type Wi vertices that are infected at time t. That is, 

Xi{t) := # {3 : w, = Wi, and iit{j) = 1} . 

Note that X(i) = {Xi{t) , X2{t) , . . . , Xm{t)) is an in-dimensional birth death process. Given the 
vertex weights w, we denote the law of X(-) started at X(0) = x by P^'^ . As with the contact 
process, we will write P^ for the law of X(-) when the associated contact process starts with 
all vertices infected (that is, when 'Yl'jLi-^ii^) — ""■)• 

The main result of this section is the following, which not only shows that the contact process 
is super-critical when A > Ac, but also gives a lower bound on the probability of a given vertex 
to be infected. 

Proposition 2.1. Let the distribution fi have finite support, and let Ac and o"(A) be defined as 
in ([6]) and ([7]). Then, if X > Xc and r] G (0, cr(A)) there exists a constant C > such that 

lim P^ ( inf Xi{t) > -— — ^pjU, Vi = 1, 2, . . . ,m ) = 1, n - a.s. 

n->oo yt<eCn l + r]Wi J 
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Proof. For any e > 0, let 

Un{e) := |x G Z'" : x, > -l^p^n , Vi = 1, 2, . . . , m 

Also, let X[ti,i2] := {^(i) : t £ [^11^2]} be the trace of the process X(f) between times ti and 
^2- The idea of the proof is that there exists 5 G {rj, o"(A)), r > 0, and C" > such that /x — a.s., 

inf P^-'' (X(r) G t/„((^) , X[0, r] C C/„(r/)) > 1 - e"^'", for ah ?x large enough. (10) 

x6C/„(<5) 

That is, we can (for n sufficiently large) uniformly bound from below the probability that 
starting from a point x G Un{S) the process X(t) a short time later is still in Un{S) and hasn't 
exited Un{f])- By dividing [0,6*-^"/^] into e*^"/^/T intervals of length r, ([10]) implies that, 
H — a.s., 

lim P^ fx[0,e^'"/2] c C/„(r/)') > lim 1 - e-^'^'/^/r = 1. 

This last statement is equivalent to the conclusion of the proposition. 

It remains only to prove (jlOll . The monotonicity of the contact process implies that we only 
need a lower bound on the inner probability in (jlOp at a single point. That is. 



inf P^'- (X(r) G UniS) , X[0, r] C C/„(r/)) = P«'-' (X(r) G [/„(5) , X[0, r] C [/„,(7?)) , 
x6C/„(<5) 

where 



1 + 5Wi 



X — l-^l) '''2> • • • ' ■''mJ' -^i •" 

For any < (5 < e let 

B„(*,.) := |x £ Z" : -i!|. < ^ < -4%. i = 1,2....™ 

[ 1 + dWi Pin 1 + eWi 

Note that Bn{S,e) C Un{S). Then, to prove (fTOj) it is enough to show that for some rj < 5 < 
e < o"(A), and constants r, C" > 0, /i — a.s., 

P^'''' (X(r) G Un{6) , X[0, r] C P„(7?, e)) > 1 - e"^'", for all n large enough. (11) 

The idea behind proving (jlip is that by choosing (5 close enough to 1] < (t(A), the drift of the 
process X(t) is increasing in all coordinates in a neighborhood of x (with high probability). 
Then, the process X(t) starting at x will with high probability increase in all coordinates by 
time r and by choosing r > small enough it will also not have exited Bn{r],£). 

As mentioned above, X(t) is an m-dimensional birth death process (given the vertex weights), 
but in order to describe the jump rates we need some additional notation. Let 



Ni,n ■■= #{j < n : Wj = Wi}, and pi,„ : = 



Ni,n 



n 

be the number and proportion, respectively, of type Wi vertices. Note that the law of large 
numbers implies that 

lim Pin =Pi, /i - a.s. 

n— >oo 

Let {ej}^]^ be the canonical basis for Z™. If X(t) = x = (xi,3;2, . . . ,Xm), then X(t) increases 
in the ith coordinate to x + ej at rate (?i^(x) and decreases in the ith coordinate to X(t) — Cj 
at rate q^ (x) , where 

g+(x) = {Ni^n - Xi) ^ Xj '-^, and gr(x) = Xj. (12) 

The following Lemma allows us to bound the jump rates in some Bnirjje). 
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Un{S) 




,X(r) 



S„(?7,e) 



Figure 1. An example of the event in (|lip . The key is to choose tj < 6 < e so 
that the drift of the process X(i) is increasing in ah coordinates inside of Bn{'i], £)■ 

Lemma 2.2. For any rj < o"(A), there exists ane € (r/,o"(A)) and positive constants {af , a~}^^ 
such that fj, — a.s., the event 

m 

n n {Q7{^)<(^7n<afn<q+{^)} 

i=lxeB„{ri,e) 

occurs for all n large enough. 

Postponing for now the proof of Lemma 12.21 we continue with the proof of (llip . Let X(t) = 
{Xi(t), . . . , Xfn{t)), where the Xi(t) are independent continuous time simple random walks that 
increase by 1 at rate afn and decrease by 1 at rate a~n. Then (j2.2p implies that (for n large 
enough) X(t) stochastically dominates X(t) while inside of Bn{ri, e). Let P^ be the law of X(i) 
started at X(0) = x. Since X(t) is a continuous time simple random walk with positive drift in 
all coordinates and jump rates proportional to n, standard large deviation estimates yield the 
following Lemma whose proof we omit. 

Lemma 2.3. Let a^ < af for i = 1,2, . . . ,m be given. Then, there exists a constant C" > 
(depending only on the a^ ) such that for any (5 > 0, 



P^ X(i)^C/„(5) <me 



-C"tn 



Vt > 0. 



On the event HilLi Clx^Bniri e) {'ii (^) < a^ n < afn < g/'(x) } , since X(t) stochastically dom- 
inates X(i) inside of Bn{rj,£) we have 

P^' (X(r) G Un{5)) < P^'"' (X(r) G C/„(5) , X[0, r] C i?„(r?, e)) + P^''^' (X[0, r] ^ i?„(r?, e)) 

(13) 
Now, the total jump rate of the process X(t) an any point x is always bounded by Cn for some 
constant C (depending on A and the Wj). Then, since the distance from x to the complement 
of Bn{rj,£) is at least cn for some constant c (depending on r],5,£ and the Wj) by choosing 
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r > sufficiently small we obtain that for some Cq > 0, 

P;j'-'(X[0,T]^i?„(7?,£))<e-^o". 
Thus, applying Lemmas 12.21 and 12.31 to (J13p we obtain that /.i — a.s., 

PZ'""' (X(r) G Un{S) , X[0,t] C Bn{r],e)) > 1 - me"^"™ - e"^"" > 1 - e-^'^ 

for all large n if C" < min{C",Co}. This completes the proof of (jlip and thus also (modulo 
the proof of Lemma 12. 2p the proof of Proposition 12.11 D 

We now return to the proof of Lemma 12.21 



Proof of Lemma \2.S[ We ultimately wish to control qfi'^) and 9j"(x) inside of some Bn{r],e), 
but fluctuations in the number of type Wi vertices force us to deal first with a slightly different 
set that depends on the actual number of type Wi vertices. For any < r] < e define 

Br,{v,e) .- |xGi^ ■ l + r^W,^ iV,,„ " I + eW,' ' ' ^'^' ' ' -"^ 

Note that the region Bn{ri, e) depends on what the actual vertex weights are, whereas the region 
Bn{rj, e) only depends on the given parameters. However, since Ni^n/n -^ Pi, IJ- — o,.s. we obtain 
that for any rj' < rj < e < e' 

Bn{j],e) C Bn{r]',e'), for all n large enough, /i — a.s. (14) 

Recall the formulas for the jump rates (?^(x) and q^{x.). Then, for any x € Bn{'r],e), 






J." 

n 



AiV,;. 



W, ™ vWf iV,„ 



^^i 



and 



1 + eWi ^^ 1 + riWj n 

_,. Wi ^. yWf 

=: 9f{r},£,n)n 



Again, since pi^n — ^ /Jj we obtain 

^lim0+(,?,.,n)=Ap.^^^^^p,^^^=:0+(ry,.), ^ - a... (15) 

and 

Jirn^ §- (?7, £, n) = pj ^ _^ ^ j^ =: Or {rj, g) , ^ - a.s. (16) 

Note that Of = {1 + A)e~ where 






i(;t 



1 + 7/^1 



10 
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Now, the definition of (t(A) implies that A(?y, ■>]) = XEfj_[wl/ {l+r]Wi)] — l > for any rj < o"(A). 
Therefore, we can choose an rj' < ij < e' < cr{X) such that A(rj',e') > 0. Fix such a ?/ and e' 
and choose af and a~ so that 

er{rj',e')<a- <at<et{r]',e'). 

Then, (^ and dH]) imply that 

9~{r]',e',n) < a~ < a^ < 9^ {rj' , e' , n) , for all n large enough, fi — a.s. (17) 

Now choose e G {ri,e'). Applying (|14p and (|17|) completes the proof of Lemma 12.21 D 

We end this section with the following Corollary. 

Corollary 2.4. Let the distribution fi have finite support. Then, for any A > Ac and r] < o"(A), 
there exists a constant c > such that 



liminf inf inf { Plimii) = 1) ^^^^ 



> 0, n — a.s. 



Proof. For A > Ac and r? < cr(A) fixed, chose c > as in the statement of Proposition 12.11 
Suppose that vertex i is a type Wj vertex. That is Wi = Wj. Then, since there are Nj^n vertices 
of type Wj , 

Therefore, 

T]Wj 



inf inf (p;j(^,(i) = l)-,^?!^^ 



inf inf 



1 



1 



K[xm 



> inf I ^£:^ 

j<m \Nj^n 

Now, Proposition 12.11 and the fact that Nj^n ~ Pjn as n — ?• 



inf X,(t) 



1 + ?7iy,- 



1 + r]Wj 



CO imply that 



lim inf ■ 



1 



-e: 



'J.fl 



i^L^.W 



t<e 



> 



rjWj 



1 + 7?W,- 



/x — a.s. 



Combining ([TH]) and ([T9|) completes the proof. 

3. Infinitely many types 



(18) 

(19) 

D 



In this section we will extend the results of the previous section to the case where fi does not 
have finite support. The first main result is a proof of part ([n]) of Theorem 11.21 

Proof of Theorem \1.2[ part ^ . For any integer tti > 1 let 

, , \xm\ 

Km[x) := A m. 

m 

Note that Km{x) < x and linim^oo f^mix) = x for all x G M. Given the sequence of vertex 

weights w = {wi}^i, let 

wM:={k^(u;.)}S,. 

The sequence vertex weights w''"^ takes on only finitely many values, and thus we may apply 
the results from the previous section. For any vertex distribution fi and m > 1, let Ac and 
(T*''"'(A) be the analogues of Ac and cr{X) for the modified vertex weights ■w'''"^ That is, 

X<r^ := {E,[K^{w,f])-\ 



and for A > Ac 

(m) 



a^^'{X) is the unique cr > that solves 1 = XEfj 



1 + aKmiwi) 



CONTACT PROCESS WITH VERTEX-DEPENDENT INFECTION RATES 11 

The monotone convergence theorem impHes that Ac \ Ac and a^"^'{X) /^ <7(A) as m — )• oo. 

If A > Ac, then there exists some m large enough so that A > Ac • For this m fixed, we may 
apply Proposition 12 . 1 1 to the modified vertex weights w^™' to obtain that there exists a C > 
such that 

lim P"(„) (r/gCn 7^ 0) = 1, ^ - a.s. (20) 

Since w'™^ < w, the monotonicity of the contact process ([8]) implies that P^{rjt 7^ 0) > 
P^^^j {r]t / 0), and thus ([20]) holds with w in place of w^™). D 

We can also extend Corollary 12.41 to the general case. This proves the lower bound needed 
for the proof of Theorem 11.31 

Corollary 3.1. Let A > Ac- Then, for any e > there exists a constant c> so that 

liminf inf inf ( PZMi) = 1) ^^^^^^] > "£, /^ - a.s. 

Proof. First note that an exercise in calculus shows that 

.^J_^.^)< --" V0<,<.. (21) 

Also, note that jr^ has derivative bounded above hy rj on x > 0, and x — Km{x) < 1/m. Thus, 
by considering separately the cases x < m and x > m 

sup — — ^^ < max <^ — , — ^ , Vr/ > 0, Vm > 1. 

Therefore, given A > Ac and e > we may choose rj < a{X) such that for all integers m large 

enough 

f_aiX}x: yumjx) \ 

sup 7-— ^-- < e. [22) 

x>o\l + a[X)x l + r]K„i{x)J 

For this choice of r/, choose m large enough so that A > Ac , V < cr(™''(A), and ()22p holds. 
Then, there exists a constant c > such that the conclusion of Corollary 12.41 holds for this r] 
and the modified vertex sequence w^™^. Then, the monotonicity property ([8]), the choice of r/ 
and m satisfying (I22p . and Corollary 12.41 imply 

liminf inf inf ( Pl(m(i) = 1) '^^^}'!'! 

> liminf inf inf ( P^,„Jr,t{i) = 1) - '^''"^'^;^ ^ - e] > -e. 



U 



4. Comparison with a multi-type branching process 



In this section we introduce a related multi-type branching process (MTBP) that we will 
use to show that the contact process dies out quickly if A < Ac. For fixed vertex weights w 
and n > 1, let Z(t) = (Zi(t), Z2(t), . . . , Zn{t)) be a MTBP with n types. Each individual dies 
at rate 1, and individuals of type i give birth to individuals of type j at rate XwiWj/n. The 
relevance of the MTBP Z(t) to the contact process r]t is given by the following lemma. 

Lemma 4.1. The MTBP Z(i) stochastically dominates the contact process rjt. That is, Z(i) 
and rjt can be coupled in such a way that tjq = Z(0) and i]t < Z(t) for all t > 0. 
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Proof. We describe the coupling informally as follows. We will divide the population of the 
branching process into two groups, one of which will be identified with the contact process 
r]t. At any time t, we will denote the individuals in groups 1 and 2 by Z^^\t) and Z(^)(t), 
respectively, so that Z(0) = Z(^)(0) + Z(2)(0). Let Z(i)(0) = r]o, so that initially all individuals 
are in group 1. The dynamics of the processes are described as follows. All individuals of type 
i in either population die at rate 1 and give birth to an individual of type j at rate XwiWj/n. 
An offspring of type j at time t is placed in group 1 if and only if the parent was in group 1 and 
Zj (t) =0 (that is, there are no individuals of type j in group 1 present already). All other 
offspring are placed in group 2. It is easy to see that Zi^^'{t) has the same distribution as the 
contact process r/t. Since Z(^)(i) > 0, the proof is complete. D 

Based on the above coupling, we will use Pw* and P^ to denote the law of Z(t) under 
the initial conditions Z(0) = e^ and Z(0) = 1, respectively. Part ^ of Theorem 11.21 follows 
immediately from Lemma 14.11 and the following Theorem. 

Theorem 4.2. Let A < Ac. Then there exists a constant C > such that 

llm P^{Z{C\ogn) = 0) = l, fi-a.s. 

Proof. Define the mean matrix M(t) = (Mij(t))"^.^-^ for the MTBP Z(t) by 

M,^,{t) = E^;;\Z,{t)]. 
Then, 

fc=i 

which in matrix form is 

-Ai = _M(t)A, 

where A = {Ai^j)"^-^ is the matrix with entries 

_ At^^iu^_ 

Solving this system of differential equations gives that 

M(t) = e^\ 

A MTBP such as Z(t) with mean matrix ]V[(t) = e^^* is called sub-critical if the largest eigen- 
value of A is negative. It is easy to see from the definition of A above that A = -■w*w„ — /, 
where w„ = {wi,'W2, ■ ■ ■ , Wn) is the row vector of weights and w* is the transpose of w„. Now, 
w*w„ is a rank one matrix with n— 1 eigenvalues at and one eigenvalue at ||w„|p = Yl^=i '^f- 
(To see this, check that w„ is an eigenvector with this as the corresponding eigenvalue.) There- 
fore, the matrix A has n — 1 eigenvalues at —1 and one eigenvalue at — 1 + -||w„|p, and the 
MTBP Z(t) is sub-critical if 

n 

A < T. TTTT = 



w^ 




For any n by n matrix Q, let ||Q|| := sup{||Qv|| : ||v|| < 1} be the spectral radius of Q. 
Then, the above analysis of the eigenvalues of A implies that 

||M(t)||=e-*(i-^'l-""'). 
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The law of large numbers and the definiton of Ac imply that 

|2 



lim 



lim — > w^ = -—, fi — a.s. 

n.->oo n ^-^ Xr- 



n— >oo n n.->oo n ^-^ Xr 

i=l 

Thus, if A < Ac there exists an e > such that 1 ||wn|P > £ for all n large enough, fj, — a.s. 

Therefore, fi — a.s., 

PZ'\Zj{t) > 1) < E^'[Zj{t)] = AUj{t) < \\M{t)\\ < e~'\ for all n large enough. 

(The second to last inequality follows from the fact that the entries of lV[(t) are non-negative.) 
Therefore, even with the initial configuration of Z(0) = 1 we have that, fi — a.s., 

n 

PZ{Z{t) / 0) < ^ PZ'\Zj{t) > 1) < n^e-'\ for ah n large enough. 

Letting t = C log n for some C > 2/e completes the proof of the Theorem. D 

We can also use the branching process approach to give an upper bound on the probability 
of a vertex to be infected after a long time. The following Proposition complements the lower 
bound in Corollary 13.11 and thus completes the proof of Theorem 11.31 

Proposition 4.3. Let A > Ac- Then, there exists a C > such that 

limsup sup sup ( P^ (7/t(i) = 1) - -— — tt^ ) < 0. 

n-!>oo t>Clogni£[n]\ l^ + 0-[A)Wi / 

Proof. The self-duality property ([9]) and Lemma |4. II imply that 

C(^*(^) = 1) = PZ'iVt / 0) < PZ'\Z{t) / 0). (23) 

To approximate Pw*(Z(i) 7^ 0) we need the following lemma. 



Lemma 4.4. For X > ( ^ E"=i % j > ^et 

a„(A) be the unique o" > that solves 1 = — > - 

n ^ 1 

Then, the extinction probability pi^n = Pi,nW satisfies 



. , - -I- awj 



/),,„ = hm PZ'' (Z(t) = 0) ^ 



Remark 4.5. If A > Ac, then an(A) and pi^n are defined for all n large enough. Moreover, 

lim a„(A) = o-(A) and lim pj „ = ^ := -— — — - — , p-a.s. (24) 

n— >oo n— >oo 1 + a[X)Wi 

Proof. An equation determining the extinction probabilities for multi-type branching processes 
is given in |AN72t Section 7.5]. In order to use this, we first need to introduce some notation. 
Let 

a. = l + E^^' m = ^ + ±^^^^s.s„ (25) 

i=i i=i 

and 

u(s) = (ni(s),'U2(s),... ,u„(s)), where ni(s) = ai(/i(s) - Sj). (26) 

Then, the extinction probability vector /5„, = {pi^n,P2,n, ■ ■ ■ , Pn,n) is the unique solution to 

u(s) = 0, < s < 1. (27) 
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Let a > and s = (si, S2, • • • , Sn), where 

1 



1 + crwi 
Recalling ()25p and ()26p we obtain that 

, n , n 

Ui{s) = 1 + --WjSj y^ WjSj -Si WiSi y^ : 

j=i i=i 

A '^ 

= (1 -Si) WiSi'y^Wj(\ - Sj) 

awi A tiij Y^ '^%' 

1 + awi n\ + awi ^-^ 1 + awj 

I \ " 2 

fjifj / A Y^ ^j 



1 + awi \ n ^-^_ 1 + aWn 

The proof is completed by noting that the term inside the parenthesis on the last line above 
equals zero when o" = ct„(A). D 

Recall that by (I23p we need an upper bound for Pw (Z(i) 7^ 0). By Lemma 14.41 we know 
that for n fixed Pw"(Z(t) / 0) decreases to 1 — pi^n as i increases. We would like to show that 
this convergence is fast enough to so that when t is large enough the error is very small. 

To this end, let Tq := inf{t > : Z(t) = 0} be the extinction time of the MTBP Z(t). 
Let PJ (•) = Pw'*(' 1^0 < 00) be the law of Z(t) started from one individual of type i and 
conditioned to eventually die out. Recall that if A > Ac then for all n large enough pi^n = 
Pw'*(ro <oo) G (0,1). Then, 

P^'*(Z(t) / 0) = Pr(Z(t) / , To < 00) + P^'*(Z(t) / , To = 00) 

< p^'XPo < oo)p:;*(z(i) / 0) + p;j'^(ro = oo) 

= pi,„P^'(Z(t) / 0) + 1 - pi,„. (28) 

Lemma 4.6. // A > Ac then there exists an e > so that, p — a.s., for all n large enough 

P^\Zj{t)^0)<e-'', yi,j<n. 

Proof. The MTBP Z(t) conditioned to go extinct can be analyzed using /i-transforms. Given 
the vertex weights w, let hn be defined by 

n 
^n(z) = JJp^;^, Z = {zi,Z2,...,Zn). 

i=\ 

Then the jump rates of the process Z(t) conditioned on extinction are given by qn{x,y) = 
g'„(x, y)/i„(y)//i„(x), where (7n(x, y) are the jump rates of the process Z(t) without conditioning 

on extinction. Therefore, we obtain that under P^ , Z(t) is a MTBP with modified birth and 
death rates. Under the measuer PJ , individuals of type i die at rate l/pi,n and give birth to 
particles of type j at rate — ^r^Pj,n- We define the mean matrix conditioned on extinction by 
M(t) = (Mi,j{t))lj=„ where 
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and E^ denotes expectations with respect to P^ . Then, as was done for the mean matrix 
M(t) above we see that 



M{t) = e^\ where A = {A,)l^^,, A, = \ x£^^ " ^''" ' ^ 



1 I XwjWj 

""" n I 

Pj,n i / J- 



Letting D^^ be the diagonal matrix with diagonal entries pi^n we obtain that 

n " ''" P" P^ \n 



:W„w:Z?^„ - D^: = D^: -i?p„w„w:Z5^„ - / 



As above, we see that the eigen values of -D^^Wny^n^pn ~ ^ ^^^ ~^ (with multiplicity n — 1) 
and -||D^^w„|p — 1 = ^ X]"=i wfpf^ — 1. Note that the defnitions of /5j^„ and CTn(A) imply that 

Therefore, since pi^n < 1 we obtain that (for n large enough) ^ XlILi '^fp'in < 1' ^^"^ thus all 
the eigenvalues of —Dp^'Wn'y^'^D^^ — I are negative. We wish, however, to show that all the 
eigenvalues of A are negative. To this end, note that D~^ is a diagonal matrix with diagonal 
entries all greater than one. We prove the following algebra lemma in the appendix. 

Lemma 4.7. Let \J he a symmetric matrix with all eigenvalues less than Aq < 0, and let D be 
a diagonal matrix with diagonal entries all greater than 1. Then, DU also has all eigenvalues 
less than Aq. 

Applying Lemma HTTl we obtain that A has all eigenvalues less than ^ Yl^=i '^Ipfn — 1 < 0. 
Recalling the definition of pi in ()24p , we obtain that 



1 " 
hm - V wfpi -l = E^ [whl] - I = E^ 



n— >oo 77, 



Wi 



(l + a(AVi)2 



1 <0, 



and therefore, p — a.s., for any < e < 1 — E^ [wl/{l + cr(A)wi)^] and all ?i large enough, 
||M(i)|| < e~^^. The proof is completed by noting that 

P^J{Z,{t)^0)<M,,{t)<\\M{t)\\. 

D 

We are now ready to finish the proof of Proposition 14.31 Recalling (123p and (j28p . Lemmas 
14.41 and 14.61 imply that, p — a.s., for all n large enough, 

P^^iVtii) = 1) < ne-^* + 1 - A,„ = ne-^' + ^ ''"i^.^f; , Vi < n. 



Therefore, recalling (|23p we obtain that 

a{X)wi 



sup sup [P^{r]t{i) = 1, 

t>Clogni6[n] V i + a[A)Wi 

f -et , 0'n{X)Wi (T{X)Wi 

< sup sup ne + 



t>C log ni£[n] 



l + an{X)wi l + a{X)wi 



<^l-Ce^ l^n(A)-<T(A)| 



V^a„(A) + ^/a{XJ 



2' 



where we used ()21|) in the last inequality. Since a„,(A) — t- o'(A) as n — t- oo, the conclusion of 
Proposition 14.31 holds for any C > 1/e. D 
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5. ASYMPTOTICS OF (7(A) 

In this section we give the proof of Proposition 11.61 The proof doesn't involve probabihty at 
all and is purely an exercise in analysis. 
Recall that cr(A) is defined for A > Ac by 



(t(A) is the unique a > that solves 1 = XE^ 



w 



Thus, it is crucial to understand the asymptotics of E^ 



1 + (TW\ 



1 + awi 
as a ^ 0+. 



Lemma 5.1. If fi{wi > x) ~ Cx (" ^^ for some C > and a > 2, then 



E„ 



wl 



and 



E,. 



1 + awi 



wi 



'c^?^a°-3 aG(2,3) 

sm(7ra) ^ V ' y 



-2Cloga 



a = 3, 



as a 



0+, 



1 + awi 



K 



/^{a-lK „_3 c,e(3,4) 

sin(7ra) ^ \ t J 

~ \ 3Ca log a q = 4 

,31 



as a 



0+. 



(29) 



(30) 



^-Ei,[wi]a a>4, 

Moreover, the conclusion in the case a > A also holds under the assumption that E^[w'f] < oo. 

Proof. We first consider the case when a < 3. Note that in this case S^[w^] = oo, and thus 
Ef^\wl/{1 + (Jt(;i)] — )• oo as (T — >• 0+. The assumption on the tail decay of fi implies that for any 
e > 0, there exists an xq = xo(e) < oo such that 



(C 



ex 



-(a-l) 



< iJ-iwi >x)< {C + e)x' 



Aa-l) 



Vx > Xq. 



Therefore, 



°° x2 



XQ 



1 + ax 



((C7-e)(a-l)x-")dx< E^ 



w\ 



1 + awi 



< 



Xn 



+ 



1 + CTXo ' Jxo 1 + <7^ 



((C + e)(a-l)x-°)dx. 



Note that t^ — and f,^" ^- — dx both are bounded as o" — t- 0"*". Therefore it is enough to show 

I+UXQ Jl 1+ax ^ 

-"-3 a E (2, 3) 
3, 



oo rj.'i-a 



I 1 + ax 



dx 



sin(7ra) 

-logo- 



a 



as a 



0+. 



a 



(31) 



When a = 3, we can explicitly evaluate this integral to be log(l + l/cr), which is asymptotic to 
— logcr as cr — 7- 0"*". To handle the case when a E (2,3) we first make the change of variables 
z = ax in the integral in (|31|) to obtain 



oo j.2-a 



dx = a' 



a-3 



oo 2:"^— a 



dz. 



Ii 1 + ax J^ 1 + z 

However, taking cr —)• 0"^ in the limit of the last integral we obtain that for any a E (2, 3) 

TT 



lim 

o--s>0+ 



oo 2— Q 



dz 



oo 2— a 



dz 



1 + z Jq 1 + z sin(7ra) 

where the last equality is a standard exercise in complex analysis. This completes the proof of 
(|31|) and thus also the proof of 
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Next, we consider the case when a > 3. In this case E^[wf] = 1/Xc < oo, and thus 
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E„ 



w 



1 + awi 



E„ 



w. 



1 



awi 



1 + awi 



1 



(tE„ 



w 



1 + awi 



(32) 



Thus, we are reduced to computing the asymptotics of E^[w'i/{l + awi)]. Note that E^[w'i/{1 + 
awi)] — )■ Ef^lwf] as (7 — )■ 0+. This completes the proof of (pO|) in the case when Ef^lwf] < oo (in 
particular when a > 4). On the other hand, if a < 4 then Efj^lwf] = oo and we must compute 
the asymptotics of how fast E^[wi/{1 + crwi)] grows as a — )- 0"'". To this end, we repeat the 
process above in the case when a E (2, 3] to show that 



E„ 



w^ 



1 + awi 



C{a-1) 



oo j.3—a 



dx 



sin(7ra) 

-3ClogCT 



a G (3,4) 
a = 4, 



as a 



0+. 



I 1 + ax 

Combining these asymptotics with ()32p completes the proof of (I30|) when a G (3, 4] . D 

With Lemma l5.ll at our disposal we are now ready to give the proof of Proposition II. 6[ 

Proof of Proposition \l.(k First, note that the definition of ct{\) and the asymptotics in Lemma 
5.11 imply that cr(A) — )• as A — )■ Ac. Then, another application of Lemma |5 . 1 1 implies that 



1 = \E„ 



wr 



'ac^^^cj(A) 



a-3 



aG (2,3) 
a = 3, 



A -Ac 
Ac 



< 



E„ 



W; 



l + a{\)wi 



~ < 



, ^ , - > as A — )• 0^ 

l + a{\)wi\ [-2CAloga(A) 

The above asymptotics are equivalent to the conclusion of Proposition [L6] in the cases a G (2, 3) 
and a = 3. Similarly, when a > 3 the defenition of cr(A) and Lemma |5 . 1 1 imply that 

' -AC;^^^cj(A)"-3 ae(3 4) 

sin(7ra) ^ ' v ' / 

-A3Ca(A)logcj(A) a = 4 asA^A+. 

^A^^[u;f]o-(A) a>4. 

The conclusions of Proposition 11.61 follow easily from the above asymptotics in the cases a G 
(3,4) and a > 4 (or E^[w''i\ < oo). In the case a = 4, the conclusion in Proposition 11.61 follows 
from the above asymptotics and the fact that 

fi5)logfi6)r^-A5 < 
where the above asymptotics are as (5 — t- 0"''. 



f{5) ~ A 



log(lM)' 



D 



Appendix A. Proof of Lemma [¥771 



Lemma 14.71 is a simple consequence (albeit not necessarily immediate) of the Spectral Theo- 
rem, but we give the proof here for completeness. 

Proof of Lemma \4.T\ As in the statement of the Lemma, let L'^ be a real nxn symmetric matrix 
with largest eigenvalue Aq < 0, and let D he an n x n diagonal matrix with diagonal entries 
Di^i > 1. Let (•,-)d be the inner product on M" defined by (x,y)£) := x*Z?~^y. Since U is 
a symmetric matrix, then the linear map x i— )■ DUx is a symmetric linear map with respect 
to the inner product {■,-)d- Therefore, by the Spectral Theorem there exists a basis {v^}^^^ 
of M" consisting of eigenvectors oi DU that are orthogonal with respect to the inner product 
{■,-)d- Moreover, by scaling the basis we may assume that ||vfc|| = 1 for all k. 

Given the basis vectors {v^} = {{vk^i,Vk^2, ■ ■ ■ ,Vk^n)}, let Ck = {'Vk,'Vk)D- Then, since the 
diagonal entries of D are all at least 1 we obtain that 



Ck 



^iD-'^k 



n 

E 

i=l 



vl ,D7} < 






1. 
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Since the v^ are eigenvectors of DU we have that DUvk = AfcVfc. Our goal is to show that 
Afc < Aq. To this end, note that 

VfcC/v/,. = {vk,DUvk)D = (vfc,AfcVfc)z) = AfcCfc. 

However, the weh known variational characterisation of the largest eigenvalue gives that 

Ao = sup{x*[/x : 11x11 = 1} > vlUvk = XkCk- 

Since c^ < 1 and Aq < this implies that A^ < Ao/cfc < Aq. □ 
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